Spectroscopy of the Schwarzschild Black Hole at Arbitrary Frequencies 
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Linear field perturbations of a black hole are described by the Green function of the wave equation 
that they obey. After Fourier decomposing the Green function, its two natural contributions are 
given by poles (quasinormal modes) and a largely unexplored branch cut in the complex-frequency 
plane. We present new analytic methods for calculating the branch cut on a Schwarzschild black hole 
for arbitrary values of the frequency. The branch cut yields a power-law tail decay for late times in 
the response of a black hole to an initial perturbation. We determine explicitly the first three orders 
in the power-law and show that the branch cut also yields a new logarithmic behaviour T~ 2e ~ 5 InT 
for late times. Before the tail sets in, the quasinormal modes dominate the black hole response. 
For electromagnetic perturbations, the quasinormal mode frequencies approach the branch cut at 
large overtone index n. We determine these frequencies up to n _5//2 and, formally, to arbitrary 
order. Highly-damped quasinormal modes are of particular interest in that they have been linked 
to quantum properties of black holes. 



The retarded Green function for linear field perturba- 
tions in black hole spacetimes is of central physical impor- 
tance in classical and quantum gravity. An understand- 
ing of the make-up of the Green function is obtained by 
performing a Fourier transform, thus yielding an integra- 
tion just above the real-frequency (w) axis. In his seminal 
paper, Leaver [1] deformed this real-w integration in the 
case of Schwarzschild spacetime into a contour on the 
complex-w plane. He thus unraveled three contributions 
making up the Green function: (1) a high-frequency arc, 
(2) a series over poles of the Green function (quasinor- 
mal modes QNMs), and (3) an integral of modes around 
a branch cut originating at u> = and extending down 
the negative imaginary axis (NIA), which we refer to as 
branch cut modes (BCMs). The three contributions dom- 
inate the black hole response to an initial perturbation 
at different time regimes. The high-frequency arc yields 
a 'direct' contribution which is expected to vanish after 
a certain finite time [21 [3] . 

The QNM contribution to the Green function domi- 
nates the black hole response during 'intermediate' times 
and it has been extensively investigated (e.g., [4j for a 
review). At 'late times' the QNM contribution decays 
exponentially, with a decay rate given by the overtone 
number n € Z + . QNMs have also triggered numerous in- 
terpretations in different contexts in classical and quan- 
tum physics, ranging from astrophysical 'ringdown' [5] to 
Hawking radiation [3 [7], the 'gauge-gravity duality' ([5] 
for Schwarzschild black holes which are asymptotically 
anti-de Sitter and [9] for asymptotically fiat ones), black 
hole area quantization [TUlfbl] and structure of space- 
time at the shortest length scales [T5]. The quantum in- 
terpretations are given in the highly-damped limit, i.e, 
for large n. The highly-damped QNM frequencies in 



Schwarzschild have been calculated up to next-to-leading 
order in [16-20 . Despite all the efforts, the leading or- 
der of the real part of the frequencies for electromagnetic 
perturbations has remained elusive (only in |21) they find 
numerical indications that it goes like n -3 / 2 ). 

The contribution from the BCMs, on the other hand, 
remains largely unexplored. The technical difficulties of 
its analysis mean that most of the studies have been 
constrained to large radial coordinate as well as small 
v = iui > along the NIA. An exception is a large-z^ 
asymptotic analysis of the BCMs in [TS] (and near the 
algebraically-special frequency in [25]) solely for gravi- 
tational perturbations. The small-;/ BCMs are known to 
give rise to a power-law tail decay at 'late' times of an ini- 
tial perturbation [TJ [21 [231 121] • In general, however, there 
is an appreciable time interval between when the QNM 
contribution becomes negligible and when the power-law 
tail starts [25] ■ The calculation of the BCMs for general 
values of the frequency (i.e., not in the asymptotically 
small nor large regimes) , to the best of our knowledge has 
only been attempted in [26, 27] where the radial functions 
were calculated off the NIA via a numerical integration 
of the radial ODE (1.2 1 followed by extrapolation to the 
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NIA, and only for the gravitational case. 

In this Letter we present the following new results: 

(1) A new analytic method for the calculation of the 
BCMs directly on the NIA and valid for any value 
of v. In particular, this method provides analytic 
access for the first time to the 'mid'-;/ regime. 

(2) A consistent expansion up to 4th order for small-f 
of the BCMs for arbitrary value of the radial coor- 
dinate. We explicitly derive a new logarithmic be- 
haviour T~ 2e ~ 5 InT at late times. 

(3) A large-f asymptotic analysis of the BCMs. It shows 
a formal divergence, which is expected to be cancelled 
out by the other contributions to the Green function. 

(4) A new asymptotic analysis for large-n of the electro- 
magnetic QNMs. The analysis is formally valid up 
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to arbitrary order in n; we explicitly calculate the 
corresponding frequencies up to n~ 5 / 2 . 
Methods in (l)-(3) provide the first full analytic account 
of the BCMs and they are valid for any spin s = 
(scalar) , 1 (electromagnetic) and 2 (gravitational) of the 
field perturbation. For the QNM calculation we focus on 
spin-1 as this is the least well understood case. We note 
that spin-1 perturbations are acquiring increasing impor- 
tance [55] [55], although it is expected that only the lowest 
overtones of the QNMs are astrophysically relevant. 

We present details in [3U] and [31]. We take units 
c = G = 2M = 1, where M is the mass of the black hole. 



I. THE GREEN FUNCTION & BRANCH CUT 

After carrying out a Fourier transform and a multipolc 
decomposition, the radial and time parts of the retarded 
Green function for linear fields on a Schwarzschild black 
hole can be written as 
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where c > 0, £ is the multipole number, r > = max(r, r'), 
r < = min(r, r') and W (w) is the Wronskian of the two 
functions fi and gi. These functions are linearly inde- 
pendent solutions of the radial ODE 
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where = r+\n(i — 1) and A = £(£+1). The solutions are 
uniquely determined when Im(u;) > by the boundary 
conditions: fg ~ e ~^r, ag r ^ _^ _ QQ an( j ^ ^ e +iur,. ag 

— > oo. The behaviour of the radial potential at infinity 
leads to a branch cut in the radial solution gi |321 133] . 



The contour of integration in Eq.( 1.1 1 can be deformed 



in the complex-w plane pQ yielding a contribution from a 
high-frequency arc, a series over the residues (the QNMs) 
and a contribution from the branch cut along the NIA: 



Gf c (r,r';t) = — / dv AG e (r,r> ; v )e~ v ' 



(1.3) 



where the BCMs are 
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with q(v) = —iAge(r,v)/gt(r,iu) where Ag e (r,v) = 
lim £ _ > .Q+ [gt{r, e — iv) — gi(r, —e — iv)] is the discontinu- 
ity of gi across the branch cut. 

We present here methods for the analytic calculation 
of the BCMs. We calculate fi using the Jaffe series, 
Eq.39 [35] . The coefficients of this series, which we denote 
by ctfc, satisfy a 3-term recurrence relation. We calculate 



gt using the series in Eq.73 [35], which is in terms of the 
confluent hypergeometric [/-function and the coefficients 
dfe. This series has seldom been used and one must be 
aware that, in order for gi to satisfy the correct bound- 
ary condition, we must set ak=o — (— 2iu) s+1 ~ 2lul , which 
itself has a branch cut. To find an expression for AGi on 
the NIA we exploit this series by combining it with the 
known behavior of the [/-function across its branch cut: 



Agi{r,v) = —, — — — x 



(r-iy r(l-2i/) 



(1.5) 



fc=0 
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r(l + s + k - 2v)T(l - s + k-2v) 



where we are taking the principal branch both for ak=a 
and for the [/-function. In order to check the convergence 
of this series, we require the behaviour for large-fc of the 
coefficients a^. Using the Birkhoff series as in App.B |34) . 



we find the leading order ak 
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calculated up to four orders higher in [3D]) as k — > oo. 
We note that this behaviour corrects Leaver's Eq.46 [32] 
in the power '— v' instead of '— 2z/. The integral test 



then shows that the series (1.5) converges for any v > 



Although convergent, the usefulness of ( 1.5 1 at small-i^ is 



limited since convergence becomes slower as v approaches 
while, for large-^, AGi grows and oscillates for fixed r 
and r' . Therefore we complement our analytic method 
with asymptotic results for small and large v. 

The small-|w| asymptotics are based on an extension 
of the MST formalism 135, 36 . We start with the ansatz 
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2-Fi(fc + fJ, + s + 1 — iu, —k — fi + s — iu; 1 — 2iu; 1 — r). 



Imposing Eq.(1.2) yields a 3-term recurrence relation for 



a k and requiring convergence as k — > ±00 yields an equa- 
tion for /x, that may readily be solved perturbatively in 
u! from starting values /z w= o = ^ an d fJ*u=o = — ^ — 1- 
Likewise for the coefficients taking ajj = 1 we obtain 



(£+l~s) 2 
2(£+l)(2£ + l) 



(£ + 1- s) 2 (£ + 2- s) 2 
~ 4(£+l)(2£+l)(2£ + 3) : 



;0J 2 + O (CJ 3 ) (1.7) 



while a_j and a^ 2 are given by the corresponding terms 
with £ — > —£—1. (Apparent possible singularities in these 
coefficients are removable.) 

The k = term in Eq.(1.6) corresponds to Page's 
Eq.A.9 [37]. To obtain higher order aymptotics we em- 
ploy the Barnes integral representation of the hyperge- 
ometric functions [38] which involves a contour in the 
complex z-plane from — ioo to ioo threading between the 
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poles of T(k+[x+s+l — iu)+z), r(— fe — n+s—iw+z) and 
r(— z). As lu — >• double poles arise at the non-negative 
integers from to max(fc + £ — s, — k — I — 1 — s), however 
we may move the contour to the right of all these ambient 
double poles picking up polynomials in r with coefficients 



readily expanded in powers of w, leaving a regular con- 
tour which admits immediate expansion in powers of v. 

By the method of MST we can also construct gg and 
hence determine q(v) and W. For compactness, we only 
give the following small-i/ expressions for the case s = 
(cases s — 1 and 2 are presented in [30 J, 
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512£ 6 + 2016^ 5 + 1616£ 4 
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where if^ is the «-th harmonic number of order r. We 
note that the In v term at second-to-leading order orig- 
inates both in q(v) and in W . In fact, both functions 
possess a In v already at next-to-leading order for small- 
i/, but they cancel each other out in q/\W\ 2 . Similarly, 
the coefficient of a potential term in g / 1 | 2 of order 
v 2 +3 (lni/) 2 is actually zero. 

Let us now investigate the branch cut contribution to 
the black hole response to an initial perturbation given 
by the field u\ c and its time derivative u\ c at t = 0: 



2Mu, 



«f c (r.,t) = (1.9) 

OO 

dr', [GfC^r'-^unri) + uf(K)d t Gf c (r,r';t)] 



OO 

We obtain the asymptotics of the response for late times 
T = t - r* - ti using Eqs.Q and (IL6|-( [L8|. We note 
the following features. The orders v 2l+z and i/ 2i + 3 in the 
BCMs AG e yield tail terms behaving like T~ 2e ~ 3 and 
rp-2i-4^ respectively. We have thus generalized Leaver's 
Eq.56 [T] to finite values of r. Furthermore, Eq.56 [T] is 
an expression containing the leading orders from uf and 
from u'f . However, the next-to- leading order from uf 
will be of the same order as the leading-order from u\ c . 
In our approach above we consistently give a series in 
small-j/, thus obtaining the correct next-to-leading order 
term for large-T in the power-law tail. Importantly, we 
also obtain the following two orders in the perturbation 
response: T~ 2i ~ 5 InT and T~ 2i ~ 5 . We note the interest- 
ing T~ 2l ~ 5 InT behaviour, which is due to the v 2l+3 \nv 
term in Eq.(1.8). To the best of our knowledge, this is 



the first time in the literature that any of the above fea- 
tures has been obtained. The logarithmic behaviour is 
not completely surprising given the calculations in [33]. 
However, one may be led to a wrong logarithmic be- 
haviour [39] if the calculations are not performed in de- 
tail. In order to exemplify our results, we give the explicit 
asymptotic behaviour in the case s = and 1 = 1 and 
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FIG. 1. (Color online). Full perturbation response ut to the 
Gaussian described above Eq.( 1.10 1 compared to the late-time 



asymptotics. Solid-red: numerical solution; dashed-black: 
Eq.(1.10l; lower curves: numerical solution minus the first 
(green), first 2 (blue) and first 4 (cyan) terms in Eq.( |1.10| . 

initial data u % [ = e - 2 (r*-xo) 2 an( j u ic _ _^ ^ — xq) u e c , 
with xq = r*(r = 5). The perturbation response due to 
the branch cut at r = 5 at late times is given by 



uf c = -2288.6 T~ 5 + 104770.5 T~ 6 + 

[86968.2 InT — 3493893.9] T~ 7 + o (T~ 7 



(1.10) 



Fig.[T]shows that these asymptotics are in excellent agree- 
ment with a numerical solution of the wave equation. 
At large-i/, we obtain the asymptotics 



2ivq(v) 



{-iy' 2 2i cos^t/) g = 

v [1 + 3cos 2 (2-7w)] ' ' ' (1.11) 
— y^iA sin(27w) 

■ W r * + sin(27rz/)e- !/r * , s = 0, 2 
fo(r, -iv) ~ { y/nX _ 
. 21/1/2 sin(27ri/) +e ' 
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These asymptotics show a divergence in Gf c when 
t < | T* | + \r'J. They also lead to a divergence in the 
perturbation response at fixed t and r for a non-compact 
Gaussian as initial data. Both types of divergences are 
expected to cancel out with the other contributions to 
the Green function. We have thus provided a complete 
account of the BCMs for all frequencies along the NIA; 
the behaviour is illustrated in Fig(2] 



II. SPIN-1 QUASINORMAL MODES 



We present here an analysis for large-n of the elec- 
tromagnetic QNMs. We may find solutions of Eq.(1.2) 
valid for fixed v e 



-1/2 



as fa = Efc=oV'i 
two independent solutions: 



~vj2 as expansion in powers of 
1.2, starting with the 
(2/v) sinhw and 



ip^'{v) = coshi>. We may express any higher order solu- 
tion in terms of the O th -order Green function as 



4 k \v) = l du s ' mh{v ~ u) 



(2w)3/2 

(2w) 1/2 [AD 2 -2D- A] 



[8L> 2 - 6D - X) 

I (fe-2)/ \ ~| 



(2.1) 



AG//M 
2 




FIG. 2. (Color online). AG i as a function of v for r» = 0.1 
and r'„ = 0.2. (a) Using Eq.(1.5l; dashed-green: s — I = 2, 
continuous-blue: a = 0, I = 1, dot-dashed-orange: s = I = 1. 
Note the interesting behaviour near the algebraically-special 
frequency |22j at v = 4 for s = 2. (b) s = 0, ^ = 1 for smal l v\ 
continuous-blue using Eq.(1.5l; dashed-red using Eq.( |1.8[ ) to 
0(u 15 ) - see 30]. (c) s = 0, £ = 1 for large ^; continuous-blue 
using Eq. ( |1.5 1 ; dashed-red using the asymptotics of Eq. ( |1.11[ ) . 



where D 



4i\e m v) 



du 



From this expression, it follows that 



_ _Jnk/2,Ik) 
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where 
arg(r) 

— iirk/i 



Oil 



— Xy/n/2. In addition, for 

tt/4, e-^-^/^iv) and 



are both real. 



It 



follows that along arg(r) = 7r/4, up to power law 
corrections, 



A z e v + B,e~ l 

1 ^ a k 
y ^- — ' i/k/2 ' 
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E 



ia x a k _ 2 



v k/2 



,fc/2 



Equating asymptotic expansions at arg(r) = 37r/4 yields 
a k e K, Im/3fc = a 2 Q!fc_2 and also serves to determine 
Re/3fc (except when k — Ap — 2 for p £ N, which do not 
contribute to the QNM condition). 



Along arg(w) = n/2 the Green functon in Eq. (2.1) is 
rapidly oscillating and we can obtain the values of a k 
directly from a stationary phase analysis 40J: 



a = 1 



Oil 



a 2 



A 2 ln2 A 
12 



a 3 



2 ' 2 
A 3 0r"(41n2-7r) 110rA 2 410rA 



48 



192 



(2.4) 
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By matching the ipi to WKB solutions along arg(r) = 
7r/4 and 37r/4 we are able to find large-i/ asymptotics for 
gi. Also, we may use the exact monodromy condition, 
ft ({r - r h )e 2n ,uj) = e 2 ™ f e {r - r h ,uj), to obtain large- 
ly asymptotics for fi. The asymptotic QNM condition 
(W = 0) in the 4th quadrant then becomes 



2{B 2 B'{ — BiB'2) ^ 
{A 1 B 2 -B 1 A 2 )B 1 f-i 

k odd 



j/fe/2' 



(2.5) 



It is straightforward to find the QNM frequencies to ar- 
bitrary order in n in terms of the a k by systematically 
solving Eq.(2.5). Explicitly, using the values in Eq.(2.4|, 
we have 



iA 2 
2n 



e - t7 r/4 7r l/2 A 3 

2n 3 / 2 



3ttA 4 
An 2 



(2.6) 



e M / 4 V^A 2 [72A 3 (tt + In 4) - 52A 2 + 41 A + 12] 
96n 5 / 2 



0(n- 3 ) 



It is remarkable that the terms in the expansion show the 
behaviour e lkn ^ 4 (R)/n k ^ 2 to all orders. In Fig.[3]we com- 
pare these asymptotics with the numerical data in |41j . 
In [31] we apply the method used to obtain Eq.( |2.6[ ) to 
the cases s = and 2 and we obtain the corresponding 
QNM frequencies up to order n -1 / 2 and have agreement 
with [HEI]. 
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